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Abstract 

We investigate the resonance type behaviour of an overdamped Brownian particle in a bistable 
potential driven by external periodic signal. It has been shown previously that the input energy 
pumped into the system by the external drive shows resonance type behaviour as function of noise 
strength. We further extend this idea to study the behaviour as function of frequency of the 
external driving force and show the occurrence of similar nonmonotonic behaviour, which can be 
ascribed as a signature of bona fide stochastic resonance. Both weak and strong driving limit 
has been explored indicating the occurrence of marginal supra-threshold stochastic resonance in a 
bistable potential system. 
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I. INTRODUCTION 



Stochastic Resonance (SR) is a nonlinear phenomena where the combined effect of noise 
and nonlinearity (bistable systems or threshold systems) leads to an enhanced response of 
a weak periodic signal with addition of noise of optimal intensity. This counter intuitive 
phenomena has been observed experimentally and in numerical and analog simulations 1|. 
However, there has been a lot of discussions in recent years regarding the validity of SR 
as bona fide stochastic resonance [2|. The archetype of SR models is represented by a 
simple symmetric bistable potential driven by a zero mean Gaussian white noise and an 
external sinusoidal bias. The response of such system has been mainly characterized by the 
response amplitude (RA) of the periodic component of the process < x(t) > or the signal to 
noise ratio (SNR). Both these response functions show typical resonance type feature with 
increasing noise amplitude D 1|. This resonance has been attributed to the matching of 
escape rate across the barrier and the frequency of the external periodic drive. An obvious 
conclusion is the occurrence of a similar peak with increasing frequency of the bias. But 
such nonmonotonic behaviour has not been observed in spite of exhaustive numerical and 
analytical studies on SNR and RA. However, ref. 3( shows the existence of peak in SNR 
with increasing frequency for a system with rectangular potential barrier and for a special 
type of driving. 

Other criteria have been proposed to study the resonance behaviour as function of fre- 
quency, thus characterizing SR as a bona fide resonance. Gammaitoni et. al. have shown 
that residence time distribution N(T) has a resonant behaviour as function of forcing fre- 
quency. Marchesoni et. al. recently showed numerically that in a Schmitt trigger N(T) 
shows a peak with frequency for both weakly and strongly driven system. The relationship 
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between SR and synchronization of 



terized by hysteresis loop 



4, 



passages from one well to another can also be charac- 



5(. This loss can be taken as a measure of SR and SR 
is shown as a bona fide resonance. Recently it is argued that input energy is also a good 
measure of SR 7|. This energy is equivalent to the work done by the external agent which 
drives the potential periodically. The input energy not only shows peaking behaviour with 
temperature (noise strength) but also takes into account only the interwell behaviour. In the 
conventional SNR both intrawell as well as interwell motion is taken into account and hence 
for small driving frequency and noise strength the motion is predominantly dominated by 
the intrawell oscillations citeiwai. Hence the peak in the input energy is a better indicator of 
matching of escape rate and the external driving frequency. This assertion has been made 
by taking into account the detailed comparison between various measures of SR. 

In this work we show that input energy not only correctly shows the matching condition 
for noise induced escape rate and the external periodic drive, but resonance type behaviour 
is also obtained as function of uj (frequency of external drive), a signature of Bonafide SR. 
Bonafide SR is a relatively new term in stochastic dynamics, by which one means the actual 
matching of the time scales T es (barrier escape time) and (period of external drive, 2ir/u). 
This is manifested as a peak in the plot of input energy with T w and D. The existence of 
SR in supra-threshold regime has also been demonstrated. We make an exhaustive study 
throughout the parameter regime of our problem and show the detailed behaviour of SR 
peak as function of D, u and A in our model. 



II. MODEL 



We consider an overdamped Brownian particle moving in a bistable potential V(x, t) = 
—x 2 /2 + x 4 /4 — Axsm(u)t) under the influence of a zero mean white Gaussian noise £(t) 
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with correlations < £(£)£(£') >= 2D5(t — t'). The Langevin equation for such system is 

x(t) = -V'(x,t)+Z(t) (1) 

and the corresponding Fokker Planck equation (FPE) is 

dP(x,t) d , Tll . ^ d , . . , . 

^Ll = .- {nx ,t)-D-)P {x ,t), (2) 

where P(x, t) is the probability density of the particle at position x at time t. The barrier 
height AV = 0.25. We consider both weak forcing Ax m < AV as well as strong forcing 
Ax m > AV limit. Depending on the parameter regime, the Brownian particle is either 
dominated by intra well or inter well oscillations. For small driving force the motion is 
mostly dominated by intrawell oscillations. At certain optimal noise strength the interwell 
motion is enhanced due to the combined effect of noise and weak periodic modulation. Such 
resonant enhancement of interwell motion is termed as Stochastic Resonance (SR). In this 
particular case the oscillatory driving force keeps the system away from equilibrium. The 
energy required to drive the system (E in ) can be calculated from Sekimoto's stochastic 
energetics formalism |6[ and its has been shown that E in depends nonmonotonically on the 
noise strength D. Since the mean current in the system is zero, hence no work is done by 
the system and all the input energy is dissipated into the bath. This suggests that input 
energy is related to dissipative loss or hysteresis loss in the system |4|, |5J]. Our effort is to 
find the signatures of SR in the input energy of the system as function of system parameters 
lj and D. 



The input energy per period(T a; ) is defined as 



7] 



E in = I < 9 q ' dt >= -Auo l < x{t) cos{ut)dt > . (3) 



t dt JtQ 

The average < .. > is done over an ensemble of particles. The asymptotic probability 



distribution P(x,t) is calculated by solving the FPE (J2J). The input energy E in can be 
rewritten in terms of asymptotic distribution as 



dx I dt x(t) cos(ut)P asy (x, t). (4) 

-oo J to 

where P asy (x, t) = P asy (x,t + TJ). This distribution is obtained after the initial transients 
have died down and asymptotically probability density assumes a unique limiting periodic 
distribution in time. The maxima of E in as function of D is taken as a signature of SR. 
We numerically solve eqn. (J2J) by method of finite difference to obtain the asymptotic 
distribution P asy (x, t). Ei n is obtained by numerically integrating eqn. (0J). Throughout this 
work all the physical quantities are in dimensionless unit 1]. 

III. RESULTS AND DISCUSSIONS 

The input energy shows a rich structure as function of A, u and D. We mainly con- 
centrate on low amplitude drive, Ax m /A < 1. In this regime input energy shows a sharp 
maxima both as function of D and to. In fig. E]we plot E in as function of D for A = 0.1 
and different values of uj. The peak in the input energy can be attributed to the syn- 
chronization of escape from the potential well and the external periodic drive as has been 



extensively discussed in previous literatures 



7[ . With increasing uj the temperature at which 



E in peaks increases as shown in fig. ^ The most favourable condition for hopping to other 
potential well is when the barrier height is minimum. When uj is increased the Brownian 
particle spends less time in the most favourable condition and hence stronger fluctuations 
are needed to cross the barrier. Barrier height also decreases with increase in amplitude of 
external drive. Hence for reasons similar to above the resonance peak shift towards lower 
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temperature as shown in the inset, where we have plotted Ei n vs D for various values of A 
as mentioned in the caption. However, our main motive is to check whether these resonance 
features are also observed when the frequency of the external drive is varied. In fig. |2] we 
plot E in as function of uo for A = 0.1 and different values of temperature. Unlike other 
response function like SNR, RA which characterize SR and has monotonic dependence on u 
, input energy shows a nonmonotonic and a peaking behaviour with u. For small amplitude 
drive and at low temperature the resonant frequency is very close to half the Kramer's rate 
for the unperturbed system ( the resonance condition being = T^/2 — > uj/tt = r&). 
For D = 0.2, the Kramer's rate (r^ = ^7=^ exp(— 0.25/ D)) is 0.0645, which is very close 
to u/n = 0.14/7T = 0.0446. With higher amplitude the resonance peak shifts to higher 
frequency for a given temperature as shown in the inset. Thus SR as a bona fide resonance 
is established. 

It is known that for A » A max , where A max is the dynamical threshold above which 
deterministic switch events take place driven by the periodic signal alone (absence of noise), 
SR is not observed. The dynamical threshold depends on both the modulation frequency 
and the wave form A(t). For the form given in eqn.( A max /A t h = 1 + f3(Q/a), where 
(3 = (2\/3/7r)gi and g\ is the smallest zero of the Airy function Ai'(—x), is order of unity 
8J]. We investigate the marginal supra-threshold regime (A > A max ) and show that SR is 
observed with increasing temperature provided A is not very large compared to A max . The 
occurrence of SR in this regime is related to noise induced stability p, |8j. In fig. El we plot 
input energy for four different valued of A at u = 1.0. Figure |U shows the corresponding 
deterministic trajectories. As shown in fig. A = 0.85 is the threshold, i.e, the particle is 
just able to cross the barrier at x — 0. Hence higher values of A is above dynamic threshold. 
For A > 0.85 the input energy shows a resonance peak with increasing temperature. This 



peak in Ei n in supra-threshold regime is shown for a very small window of A. For the specific 
parameter values as above, the peak vanishes for A > 0.91. This peak can be ascribed to 



8]. However, in supra-threshold 
9 1 for all values of A. This 



resonant trapping as discussed in previous literature 
regime E in clearly exhibits peaking behaviour as function of u 
is more akin to conventional resonance (absence of barrier). 

In conclusion we have calculated that the input energy pumped into the system by an 
external drive using the method of stochastic energetics. This input energy is shown to be 
a good quantitative measure of SR. Moreover SR is shown to be a bona fide resonance. 
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FIGURE CAPTIONS 



Fig 1. E in vs D for A = 0.1. The inset shows the plot of E in vs D for 
u = 0.1, 0.2, 0.3 from top to bottom at A = 0.2. 

Fig. 2. Ei n vs u for A = 0.1. Inset is the plot of E in vs u at A = 0.2 for 
D = 0.15, 0.2, 0.25 from top to bottom curve respectively. 

Fig. 3. E in vs D at u = 1.0 for large amplitude A = 0.85, 0.87, 0.90 
and 0.92 respectively. 

Fig. 4. The determistic trajectories (D = 0) for amplitudes 
A = 0.85, 0.90 and 0.92. All these trajectories cross the barrier at x = 0. 



9 




0.1 



0.25 



0.05 



0.15 d 



0.25 



FIG. 1: 



1.40 - 

1.15 

0.90 



0.65 [ / 

/ 




D = 0.15 
D = 0.20 
D = 0.25 



0.40 



0.05 



0.15 



0.25 



0.35 



CO 



FIG. 2: 



10 



2.1 



2.0 



1.9 




FIG. 3: 



0.5 



0.25 



0.25 



-0.5 



0.75 



1.25 



4985 



4990 



4995 



5000 







/ \ 
/ \ 
/ \ 




/ 
/ 


\ 
i 




,-\ 
/ \ 
/ \ 
/ i 






if \ \ 




-h 


\\ 




/' N \ : 
>/ > i 




1 
II 








»'. 


-J 


! v , 


1 

1 
1 
1 

— 1 






1 










1 
1 
1 
t 




n 


1 

1 
1 








— fi- 

i 


1 

1 t 

\\ 




I 


W 




1 1 i 


/ 




35 




1 






.\\ 


/ 


\ 



4985 



4990 
Time t 



4995 



5000 



FIG. 4: 



11 



